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We describe a measurement of the top quark mass using events with two cliarged leptons collected 
by the CDF II detector from pp collisions with = 1.96 TeV at the Fermilab Tevatron. The likeli- 
hood in top quark mass is calculated for each event by convoluting the leading order matrix element 
describing — > ^ hlvibl'i/gi with detector resolution functions. The presence of background 
events in the data sample is modeled using similar calculations involving the matrix elements for 
major background processes. In a data sample with integrated luminosity of 340 pb~^, we observe 
33 candidate events and measure Mtop = 165.2 ± 6.1(stat.) ± 3.4(syst.) GeV/c^. This measurement 
represents the first application of this method to events with two charged leptons and is the most 
precise single measurement of the top quark mass in this channel. 

PACS numbers: 14.65.Ha, 13.85.Ni, 13.85.Qk, 12.15.Ff 



I. INTRODUCTION 



The standard model accommodates quark masses 
through Yukawa couphngs to the Higgs boson, but does 
not predict the size of these couphngs and contains no 
explanation for the observed quark masses. A striking 
feature is the large mass of the top quark, the heaviest 
of the observed fundamental particles. Its large mass 



suggests that it maj^lay a unique role in electroweak 
symmetry breaking P, Q ■ Precise measurements of the 
mass of the top quark constrain the mass of the yet unob- 
served Higgs boson through radiative corrections Q , and 
can restrict possible extensions to the standard model Q . 

In collisions of protons and anti-protons at ^/s = 1.96 
TeV, top quark pairs are produced primarily through the 
annihilation of quarks and anti-quarks. In the standard 
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model, the top quark decays to a 6 quark and a W boson 
nearly 100% of the time; the W boson decays to a pair 
of quarks or a charged lepton and neutrino. Quarks frag- 
ment and hadronize and are reconstructed as jets (clus- 
ters of particles). The dilepton channel, consisting of 
decays ti biviM' , has a small branching fraction, 
but measurements of the mass in this channel have the 
advantage that they are less reliant on the calibration of 
the jet energy scale than channels with hadronic W boson 
decay. A top quark mass measurement in this channel is 
an important verification that the observed top quark 
candidates are consistent with standard model produc- 
tion and decay. A discrepancy from measurements in 
other channels could indicate the presence of physics be- 
yond the standard model that makes contributions to the 
dilepton sample j^. 

Reconstruction of the top quark mass in the dilepton 
channel is particularly challenging due to the two un- 
detected neutrinos from the W boson decays. Previous 
measurements 0, Q in this channel using Tevatron run 
I data calculated a mass in each event by making sev- 
eral kinematic assumptions and integrating over the re- 
maining unmeasured quantities; the distribution of event 
masses was then compared to simulation at varying top 
quark masses. This article describes in detail the first 
application to the dilepton channel of a technique pi- 
oneered for analysis of single lepton tt Mvjjqq' de- 
cays H, 1^ 0, Q • This technique convolutes the ma- 
trix element for tt decays with detector resolution func- 
tions and integrates over unmeasured quantities to con- 
struct per-event likelihoods in top quark mass. We relax 
many of the kinematic assumptions of previous meth- 
ods and integrate over six unmeasured quantities. The 
event likelihoods are directly multiplied to obtain the 
joint likelihood from which Mt is extracted. This weights 
events according to the relative amount of information 
they carry. The data used in this measurement corre- 
spond to an integrated luminosity of 340 pb~^ collected 
between March 2002 and August 2004 by the CDF II 
detector in colhsions at ^/s = 1.96 TeV at the Fermilab 
Tevatron; this measurement and a brief description were 
first reported in Ref . ■ 

Sections II and III describe the CDF detector and the 
selection of the data sample. Section IV gives an overview 
of the analysis method. Sections V, VI and VII describe 
in detail the major pieces of the likelihood calculation. 
Section VIII describes the reconstruction and calibration 
of the top quark mass. Section IX covers the systematic 
uncertainties and Section X presents the measurement. 



II. THE CDF DETECTOR 

The CDF II detector is an azimuthally and forward- 
backward symmetric detector designed to study pp col- 
lisions at the Fermilab Tevatron. The CDF coordinate 
system is right-handed, with the z-axis pointing along a 
tangent to the Tevatron ring along the proton direction. 



The remaining rectangular coordinates x and y are de- 
fined pointing outward and upward from the Tevatron 
ring respectively, and the azimuthal angle 4> is measured 
relative to the x axis in the xy-plane. Transverse quan- 
tities such as transverse momentum, pT, and transverse 
energy, Et^ are projections onto this plane. The polar 
angle 6 is measured from the proton direction and is typ- 
ically expressed as pseudorapidity rj = — ln(tan |). Sub- 
detectors which are particularly relevant to this analysis 
are described below. A more complete description of the 
CDF II detector can be found elsewhere [T^ . 

The CDF tracking system consists of an inner sili- 
con microstrip detector and a large outer open-cell drift 
chamber. These subsystems are immersed in a super- 
conducting solenoid producing a 1.4 T magnetic field 
parallel to the p and p beams. The silicon detector, 
which provides high-resolution position measurements of 
charged particles close to the interaction region, consists 
of three subdetectors. The innermost detector. Layer 00 
(LOO) ^3 J is a single-sided layer of silicon wafers mounted 
directly on the beampipe at a radius of 1.6 cm. The 
SVXII il6j detector is a five layer, double-sided silicon 
detector that covers the radial region between 2.5 cm 
and 10.6 cm. The Intermediate Silicon Layers (ISL) 
comprise one or two additional layers of double-sided sil- 
icon, depending on the polar angle, at radii from 20 cm 
to 28 cm. The Central Outer Tracker (COT) [ll|, a large 
open-cell drift chamber, is positioned outside the silicon 
detector from radii of 0.43 m to 1.32 m. The COT con- 
tains 8 superlayers (alternating between axial and ±2° 
stereo angle) each containing 12 wire layers for a total of 
96 layers. In combination the Silicon and COT detectors 
provide excellent tracking up to \ri\ < 1.1 with decreasing 
coverage to \ri\ < 2.0. 

Sampling calorimeters segmented in r] and (j) surround- 
ing the tracking system measure particle energies. In the 
central region of |?7| < 1.1, the calorimeter is divided 
into projective towers subtending 15° in (j) and 0.1 in 
7]. The central electromagnetic calorimeter (CEM) 
constitutes the front of the wedges in the central region. 
The CEM consists of alternating layers of lead and scin- 
tillator, amounting to 18 radiation lengths of material. 
Embedded in the CEM is the shower maximum detector 
(CES). The CES provides position measurements of the 
electromagnetic showers at a depth of 5 radiation lengths 
and is used in electron identification. Behind the CEM is 
the central hadronic calorimeter (CHA) which pro- 
vides energy measurements of hadronic jets. The CHA 
consists of 4.7 interaction lengths of alternating steel and 
scintillator. In addition to the central calorimeters, end 
plug calorimeters cover 1.1 < \r]\ < 3.6. The plug electro- 
magnetic calorimeter (PEM) j21| consists of alternating 
lead absorber and scintillating tile readout with wave- 
length shifting fibers; the total thickness is 23.2 radiation 
lengths of material. A plug shower maximum detector 
(PES) 01 provides position measurement of electron and 
photon showers. The plug hadronic calorimeter (PHA) 
has alternating layers of iron and scintillating tile for a 
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total of 6.8 interaction lengths. 

The muon detection system consists of three sand- 
wiched drift tube layers, each utilizing single wire 
drift cells four layers deep. Directly behind the cen- 
tral hadronic calorimeter is the central muon detec- 
tor (CMU) which can detect muons with pT > 
1.4 GeV/c in the region of \ri\ < 0.6. Additional muon 
coverage in this region is provided by the central muon 
upgrade (CMP) which is separated from the CMU by 
60 cm of steel. The CMP detects muons with px > 
2.0 GeV/c. The central muon extension (CMX) provides 
further coverage in the region of 0.6 < jryj < 1.0. 

CDF's three level trigger system reduces the event rate 
from 1.7 MHz to « 80 Hz. The first two levels are hard- 
ware triggers that partially reconstruct events using in- 
formation from individual subdetectors while the third 
level is a software trigger that performs event reconstruc- 
tion. 



III. DATA SAMPLE 

We select ti MviM'Vg' decays with a high-p^ lepton 
trigger and the requirement that candidates have (i) two 
leptons each with px > 20 GeV/c, (ii) significant miss- 
ing energy transverse to the beam direction {^t) [23 |. 
and (iii) two jets each with Et > 15 GeV. The selec- 
tion was designed for a cross-section measurement and is 
described as "DIL" in Ref. [2^ . A description of the trig- 
ger requirements and selection used to obtain this dataset 
follows. 



A. Trigger 

The trigger requires at least one high-p^ lepton. For 
central electron candidates, the first two trigger levels re- 
quire an electromagnetic calorimeter cluster with a con- 
firming COT track and without a large hadronic energy 
deposit. The third level trigger requires an electron can- 
didate with Et > 18 GeV. Events with electron can- 
didates in the plug (I77I > 1.2) are required to have 
electron Et > 20 GeV and missing transverse energy 
^T > 15 GeV. For muon candidates, the first two trigger 
levels require hits in the muon chambers and a confirm- 
ing COT track. The third level trigger requires a muon 
stub with a matching track of pt > 18 GeV/c. 

B. Leptons 

Final lepton requirements are tighter than those made 
in the last stage of the trigger. Electron candidates are 
required to have an electromagnetic calorimeter cluster 
with Et > 20 GeV and \r]\ < 2.0. Muon candidates are 
required to have a track with pT > 20 GeV/c, which lim- 
its them to the rj coverage of the COT. At least one of the 
leptons is required to be isolated. Electrons and muons 



are isolated if the energy deposited in the calorimeter in 
a cone AR = ^(A?y)2 + (A0)2 = 0.4 around the lep- 
ton is at most 10% of the lepton transverse energy. In 
addition, electron candidates are required to have a well- 
measured track pointing at an energy deposition in the 
calorimeter. For electron candidates with \ri\ >1.2, this 
track association uses a calorimeter-seeded silicon track- 
ing algorithm . Muon candidates are required to have 
a well-measured track linked to hits in the muon cham- 
bers and energy deposition in the calorimeter consistent 
with that expected for muons. If the event contains two 
muons, only one is required to have hits in muon cham- 
bers used in the trigger decision. The other muon may 
have hits in chambers not used for the trigger decision 
if there is a matching COT track, or no hits in muon 
chambers if the COT track points in regions where there 
is no muon chamber coverage. 

C. Jets 

A jet is defined cluster of energy surrounding a 
calorimeter tower with Et > 1 GeV, grouped within a 
cone of AR = ^/{Atj)'^ + (A(/))2 = 0.4 by the jetclu 
algorithm [23|. Events are required to have at least two 
jets with I77I < 2.5 and Et > 15 GeV after the corrections 
described below are applied. 

The total jet Et is corrected for non-uniformities in 
the response of the calorimeter as a function of j], for 
effects from multiple pp collisions, and for the hadronic 
jet energy scale of the calorimeter The two highest 
Et jets for each event are assumed to stem from the b 
quarks; this assumption is true for ~70% of simulated tt 
events. The momentum components of each b quark are 
then calculated from the measured jet Et and angle by 
assuming a b quark mass of 4.7 GeV/c^ No explicit 
identification of b jets is used. 



D. Final Selection Cuts 

After lepton and jet identification, further require- 
ments are made to reduce the expected level of back- 
ground in the sample. Events are required to have miss- 
ing transverse energy of ^t > 25 GeV. ^t is corrected 
for the presence of isolated high-py muons by subtracting 
the momentum lost by the muons in the calorimeter and 
adding the muon pT to the vector sum. In events with 
^T < 50 GeV, the direction of the ^t vector is required 
to be separated by at least 20° in (p from any lepton or 
jet in the event. This reduces the background from Drell- 
Yan production of r pairs as well as the number of events 
in which mismcasured jet or lepton energy contributes a 
large fraction of the ^t ■ 

To reduce the number of Z/^* — > ee, fifj. events 
in which mismeasured jet energy leads to significant 
amounts of measured missing tranverse energy, ee and 
fifi events with dilepton invariant mass between 76 and 
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106 GeV/c^ are required to have their v ector point 
away from any energetic jets in the event p30l |. 

To further suppress background, events are required to 
have Ht; defined as the scalar sum of lepton Et, the Et 
of the two leading jets, and ^t, greater than 200 GeV. 

Events which are likely to be due to cosmic rays are 
removed by requiring a coincidence of the muon arrival 
times to the calorimeter. Electrons resulting from photon 
conversion to e^e~ pairs are also removed. Conversions 
are identified by pairing the electron track to a track 
of opposite sign and requiring that the two tracks are 
consistent with originating from a common vertex and 
being parallel at that vertex. Events with three leptons 
are removed as well as events in which the leptons have 
the same sign. 



E. Sample Composition 



TABLE I: Expected numbers of signal and background events 
for a data sample of J Cdt = 340pb~^. The signal cross sec- 
tion is obtained from Isill . The total expected background is 
the sum of the indented background contributions. 



Source 


Events 


tt (Mt = 165 GeV/c^ a = 9.1pb) 
tt (Mt = 175 GeV/c^ o = 6.7pb) 
tt (Mt = 185 GeV/c^ cr = 4.9pb) 


21.7 ± 1.4 

17.2 ± 1.4 

13.3 ± 1.4 


Total Expected Background Rate 


10.5 ±1.9 


WW 


1.2 ±0.2 


WZ 


0.4 ± 0.1 


Z/'y* ee, fj,^ 
Z/-y* — > TT 
Lepton misID 


4.7± 1.2 
0.8 ±0.2 
3.5 ± 1.4 



Total Expected Rate (Mt = 165 GeV/c^) 32.2 ± 2.3 

Total Expected Rate (Mt = 175 GeV/c^) 27.7 ± 2.3 

Total Expected Rate (Mt = 185 GeV/c^) 23.8 ± 2.3 

Run II Data (/ Cdt = 340pb"^) 33 



Tablc^lists the number of expected background events 
of each type and the number of ti signal events expected 
at various top quark masses jsj for the data sample used 
in this measurement. The signal estimate includes tt 
events in which a W boson decays to a r when the r 
decays to an e or a ^. Studies in Monte Carlo simula- 
tions show that 14% of the accepted signal events have 
at least one W boson decaying to a r. 

The largest source of expected background is Z/-f* 
ee, nil events with associated jets. The expected contri- 
bution to the sample from this background is estimated 
using a combination of Z boson data and PYTHIA 
Monte Carlo events. The second largest source of ex- 
pected backgrounds is events in which one or more jets 
are misidcntified as leptons. The probability for a jet to 
be misidcntified as a lepton is very small, so the majority 
of these events have one real lepton and one jet misidcnti- 
fied as a lepton. The expected contribution to the sample 
from this background is estimated using W £i/+}ets 
data. There arc several smaller sources of expected back- 
ground: WW,WZ and Z/j* tt with leptonic r de- 
cays. The expected contribution from these processes is 
estimated using ALPGEN [s^ and PYTHIA Monte Carlo 
events. Other processes make negligible contributions. 
A detailed description of the method of background esti- 
mation used can be found in Ref. 1251. 



IV. ANALYSIS OVERVIEW 

The probability density for ti decays is expressed as 
Ps(x.\Mt), where Mt is the top quark pole mass and x 
contains the lepton and jet momentum measurements. 
We calculate Ps(x.\AIt) using the theoretical description 
of the tt production and decay process expressed with 
respect to x, 



where ^ is the differential cross section evaluated with 
respect to event measurements contained in x. 

To evaluate the differential cross section ^ , we convo- 
lute the leading order matrix-element M for qq ~* tt ^ 
blvibl'vi' with detector resolution functions and integrate 
over unmeasured quantities. 

The matrix element depends on the momenta of the 
incoming partons (qi and 52)1 and of the outgoing two b 
quarks [pi andp2), two leptons (£1 and £2), and two neu- 
trinos [vi and ^2). Observed quantities consist of jets ji 
and j2, measured leptons Li and L2, and two components 
of missing transverse energy. To express the differential 
cross section with respect to these observed quantities x, 
transfer functions are introduced to connect the quan- 
tities which correspond to external legs of the matrix 
element {qi,q2,Pi,P2,hJ2,i'i,'^2) to the observed quan- 
tities (ji , i2, Li, L2). Quantities which are well- measured 
by the detector, lepton momenta and jet angles, are de- 
scribed by delta functions which directly reduce the num- 
ber of unknown parton-level quantities. Integrations are 
performed over quantities which are not directly mea- 
sured, i.e. quark and neutrino energies. While quark en- 
ergies arc not directly measured, they can be estimated 
from the observed energies of the corresponding jets. The 
transfer function between quark and jet energies param- 
eterizes this relationship, and is expressed as W(Ep, Ej), 
the probability of measuring jet energy Ej given parton 
energy Ep. 

We make the following assumptions regarding the 
transfer between parton-level quantities and the observ- 
ablcs: 

• Leptons are measured perfectly. We express the 
lepton transfer functions as a three-dimensional S- 
function, 

5^£i-Li)S\£2-L2). 



(1) 



• Jet angles are measured perfectly, and jet energy 
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can be described as a parametric function of par- 
ton energy. We express the b quark to jet transfer 
function as 

• The two most energetic jets in the event are due to 
the fragmentation and hadronization of the two b 
quarks from top quark decay. All other jets in the 
event, if present, are disregarded. 

• Incoming partons arc massless and have no trans- 
verse momentum. 

• Masses of the final state leptons are zero, masses of 
the b quarks are set to 4.7 GeV/c^. 

The probability density in x for qq ti b£~h'ibi''^h''g 
for a fixed Mt can be written as 

PMMt) = J d^Y.\^tliQ^.P^■,Mt)\^ (2) 

a,b 

xl[W{Ep,E,)ff^,{q,)f^M- 

jets 

In this expression, the integral is over the phase space 
(i$ for qq ^ ti blvibl'vii, the sum runs over the flavors 
a, 6 of the incoming partons, and /p^p are parton distri- 
bution functions for flavor a. Constraints such as con- 
servation of momentum which appear as delta functions 
and modify the integration are here implicitly included 
in the phase-space integration and are discussed in detail 
in the following sections. The term, \/a{Mt), in front of 
the integral ensures the normalization condition for the 
probability. 



B. Calibration 

The probability P(x) is an approximation of the true 
probability and balances precision with computational 
tractability. To account for the approximations made 
in the construction of -P(x), we derive a calibration in 
fully realistic Monte Carlo events. This strategy accounts 
for the differences between the model used in construct- 
ing P(x) and the fully realistic simulation. Uncertain- 
ties in the model used to produce the simulated events 
contribute to systematic errors on the final measurement 
described in Section VIII. 



V. TRANSFER FUNCTIONS 

The transfer function between quark and jet energies, 
W{Ep, Ej), expresses the probability of measuring jet en- 
ergy Ej from a given parton with energy Ep such that 



n{Ej,Ep)dEjdEp = n{Ep)dEpW{Ep, Ej] 



(5) 



where n{Ej, Ep)dEjdEp is the number of events with jet 
energy between Ej and Ej + dEj and parton energy be- 
tween Ep and Ep + dEp, and n{Ep) is the number of 
partons with energy between Ep and Ep + dEp. 

We parameterize the distribution of measured jet ener- 
gies, Ej, as a function of the quark energies, Ep, and the 
difference between the parton energy and the jet energy, 
W{S = Ep ~ Ej) The parameterization is a sum 

of two Gaussians to account for both the peak of the 5 
distribution and its tails, 



dx Ps{-x\Mt 



(3) 



1 



where the integration is performed over all accepted x to 
account for mass-dependent effects of the selection. 



A. Signal and Background Processes 

We calculate the probability for the dominant back- 
ground processes, Pbg(x) and form the generalized per- 
event probability density in x, 

P{yi\Mt) = Ps{^\Mt)p.,{Mt)+Pbg,{lL)pbg,+Pbg.{yi)Pbg.+- ■ 

(4) 

as a weighted sum of the probabilities for each process, 
where the weights Ps{Mt) and pbg- are determined from 
the expected fractions of signal and background events 
(see Table IJ). We evaluate probabilities for the three 
largest expected backgrounds: Z/j* boson production 
with 2 associated jets, W boson production with 3 asso- 
ciated jets in which one jet is incorrectly identified as a 
lepton, and W boson pair production with 2 associated 
jets. 



W{S) = 

V2'k{p2 +P3P5 

where each pi depends linearly on Ep: 
Pi^ ai + biEp. 



-(<S-P4)^ 
2p? 



, (6) 



(7) 



The parameters pi arc extracted with an unbinned 
maximum likelihood fit over N jets in a sample of simu- 
lated HERWIG [131 tt events with Mt = 178 GeV/c^ which 
pass the event selection and contain jets whose axis is 
contained in a cone of Ai? = 0.4 surrounding b quarks. 
Jets which arise from initial or final state radiation are 
excluded. The log likelihood is expressed as a sum over 
jets: 



-InL = - 



TV 

E 



N 



ln,i(i?pJ-^lnVK(i?p,,i?,J. (8) 



k=l 



The first term does not depend on the parameters pi 
and can be dropped from the minimization. We extract 
the parameters shown in Table UTI 
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TABLE II: Parameters for W(Ep,Ej) extracted using jets 
matched in angle to b quarks (see text) , from herwig tt Monte 
Carlo. 



Pi Oi bi 

Pi 1.90 ± 0.62 GeV 0.023 ± 0.008 

P2 2.83 ± 0.54 GeV 0.075 ± 0.005 

P3 0.70 ± 0.08 0.000 ± 0.001 GeV"^ 

P4, -1.79 ±0.79 GeV -0.187 ±0.012 

P5 8.04 ±0.67 GeV 0.095 ± 0.008 



process. We confirm this by using W{Ep, Ej) parame- 
terized from Mt = 178 GeV/c^ events to calculate the 
jet energy distribution resulting from b quarks in Monte 
Carlo top quark decays of varied top quark masses. Fig- 
ure 13 shows that the jet transfer function derived using 
partons from Mt = 178 GeV/c^ top quark decays satis- 
factorily describes jet energies from top quark decays of 
Mt ranging from 150 GcV/c^ through 200 GcV/c^. The 
performance of the transfer functions in fully realistic 
simulated events is included in the measurement calibra- 
tions discussed below. 




"0 20 40 60 80 100 120 140 160 180 200 220 240 
E, (GeV) 




20 40 60 80 100 120 140 160 180 200 220 240 
E, (GeV) 



FIG. 1: Top, difference of Ep — Ej between parton and jet en- 
ergy for reconstructed jets matched in angle to b quarks (see 
text). Center, distribution of Ej of jet energy. Bottom, input 
distribution of parton energy Ep. Histograms are simulated 
events for Mt — 178 GeV/c^; curves in the upper two his- 
tograms show the distributions calculated using W{Ep,Ej). 



To test the jet transfer function, we calculate the dis- 
tribution of jet energies which result from simulated par- 
tons of known energy. The calculation for jet energies 
resulting from partons with energy Ei < Ep < E2 is the 



integral of n{Ej, Ep)dEp: 



El 



n{Ep)dEpW{Ep,Ej). 



(9) 



VI. SIGNAL LIKELIHOOD 

The probability density for qq it s- blvihi'vii decays 
is constructed as the differential cross section, dcr, with 
respect to the measured event quantities, x. The total 
cross section a is written as 



E 

a.b 



{2nr\M\ 



2 



^ /pDI 



Will ■ q2Y 



./PDF 



Ehe an 



(10) 



d^(idqidq2, 



where the sum runs over incoming parton flavors, M is 
the matrix element for the process, gi.2 and toi,2 refer to 
the momenta and mass of the incoming partons, fpoF 
are the parton distribution functions for flavor a, and 
the integration is over the phase space for the six final 
state particles as well as the longitudinal momenta of the 
incoming particles. 

The matrix element Ha has the form 



\M\' 



FF ((2 - P^sD - X. 



(11) 



where [3 is the top quark velocity in the qq rest frame, Xsc 
contains terms describing spin correlations between the 
top quarks, Qs is the strong coupling constant {gl/Air = 
Qfs), Sqt is the sine of the angle between the incoming par- 
ton and the top quark, and F and F are the propagators 
for the top and the anti-top quark respectively. We drop 
the spin correlation term X^c as it is negligible. The top 
quark propagator and decay terms are given by 



9n 



"Iv 



(m 



(A'/tFt 



(12) 



™?(1-C|J+™| Jl+Cf-fe) 



M2,)2 + {MwVwf 



The calculation for the jet energy distribution and the 
difference in jet and parton energies resulting from all 

partons ( < iJp < 1 TeV ) in a simulated sample of where nit is the invariant mass of the t quark decay prod- 
ti with Mt = 178 GeV/c^ is shown in Figure [T] Similar 
tests using slices of parton energies are shown in Figurc|21 
The jet transfer function models the detector response 
to partons and should be independent of the production 



nets and Cy is the cosine of the angle between particles i 
and j in the W boson rest frame. The Mj, F^, M^^, F^^ 
arc the pole masses and widths of the top quark and W 
boson, and g^, is the weak coupling constant. The top 
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FIG. 2: Comparison of simulated Ej with calculations from W{Ep, Ej) from six ranges of Ep. Histograms are simulated events; 
curves show the calculated distributions using W{Ep, Ej). 
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FIG. 3: Comparison of simulated Ej with calculations from W{Ep, Ej) from distributions Ep in simulated samples with Alt ~ 
150, 160, 170, 180, 190, and 200 GeV/c^ Histo grams are simulated events for Mt = 178 GeV/c^; lines show the calculated 
distributions using W{Ep,Ej) derived using partons from Mt = 178 GeV/c^. 



quark width, Ff, is a function of A^t,Mw and Tw 
as described by the standard model. F is given by the 
same expression as Eq. 1121 replacing the terms for t and 
its decay products with i and its decay products. 



While approximately 15% of tt pairs in pp collisions 
at = 1.96 TeV are produced in gluon-gluon fusion 
(gg —> ti) |37j . our studies have shown that this term can 
be excluded from the matrix element with very little loss 
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of sensitivity to the measurement. A small systematic 
uncertainty is derived from theoretical uncertainty in the 
relative gluon fraction; see Section IX. C. 

To evaluate the differential cross section with respect 
to observed quantities, we introduce conditional 

probability terms that relate the observed quantities to 
the parton level variables and subsequently integrate 
over unconstrained parton-levcl quantities, as described 
above. 

The differential cross section is given by 



da ( 1 



/ ^ |M(gi,g2 ^ 4,i^i,Pi,4,J^2,P2) 



fe=1.2 

X ^{qxi)S{qyJdq.^, dqy,6{q^^) d{qy,J dq^^ dqy^ 
^ (27r)VpV(g|yi^t,ea.) fUilijE^eMl dqj 
A J [qi ■ (72)^ ~ "li^, 



d^i^i d^V2 



2Ey-^ 2E^2 



(13) 

The sum is over possible incoming parton fla- 
vors. The term l/4y^ {qi ■ 52)^ — '^qi"^q2 reduces to 

/flux (921, 922) = l/2'Z2i'722 with the assumption that the 
parton mass is small in comparison to the longitudinal 
momentum. 

The measured missing transverse energy is a sum of 
real missing transverse energy due to neutrinos, jet en- 
ergy mismeasurement and unclustered energy from soft 
recoil which are not included in the reconstructed ob- 
jects. The contributions from neutrinos and jet energy 
mismeasurement are accounted for by the matrix element 
and the transfer functions, respectively. Therefore, differ- 
ences between the measured and solved (as below) miss- 
ing transverse energy, J^t, are attributed to unclustered 
energy in the calorimeter. The /utf factor describes the 
prior probability of observing unclustered energy in the 
event; it is parametrized as a Gaussian in each of the x 
and y axes with no correlation. A width of 12 GeV/c^ is 
extracted from a fit to simulated samples. 



A. Phase Space Transformation and Integration 

We integrate over the lepton momenta, initial parton 
momenta, intermediate top quark and W momenta, an- 
gular components of the b partons and the six compo- 
nents of neutrino momenta. 



In order to efficiently integrate over the parton-level 
variables, we perform a transformation which splits the 
original phase space into subspaces and introduces the 
equivalent number of extra variables and integrations. 
We introduce invariant masses that correspond to inter- 
mediate t and t quarks and W bosons. Each additional 
integration over an invariant mass of the intermediate 
particle has a corresponding (5-function in squared invari- 
ant mass, and each intermediate particle four-momentum 
has corresponding (5^-function for the momentum conser- 
vation at the intermediate vertex. 

The expression of the integrand is written in terms of 
the momenta of the final state particles, Mvfhl'v'^. Inte- 
gration over the t quark and W boson invariant masses 
(mt . and mv^i ) requires expressing the neutrino momenta 
in terms of these invariant masses. These two sets of vari- 
ables are related by a system of six coupled quadratic 
equations written in terms of the final-state momenta 
and the W boson momenta {Wi) and derived from ex- 
pressions in the 5 functions. Note that we explicity as- 
sume that the ti system has no transverse momentum 
and therefore do not use the measured missing tranverse 
energy to derive the solutions for neutrino energies. Solu- 
tions in which the measured missing transverse energy is 
significantly different from the sum of the neutrino tran- 
verse energies arc dewcightcd by the unclustered energy 
transfer function, /utf. 



ml = (pi + W^f 
ml = {p2 + W2f 



1^1 1 



V2] 



(14) 



{pi+ll+Vl +P2+l2+V2)x = 
(pi + £1 + Z^i + P2 + ^2 + y2)y = 



We rewrite these equations as a single fourth-order 
polynomial and find the solutions numerically using the 
Sturm Sequence approach |38l |. The transformation be- 
tween the phase space for neutrino momenta and invari- 
ant masses is not one-to-one due to the non-linearity 
of the relations. Multiple neutrino solutions may exist 
for specific invariant masses; such solutions are there- 
fore summed. Other invariant masses may have no cor- 
responding region of neutrino phase space and therefore 
no solutions and no contribution to the total probability. 

Finally, the transformation of variables requires the 
inclusion of a Jacobian term, J. The final form of the 
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expression is 



da 



/ 1 



(2-)^ / E 



J- 



a.b 



\pi 1 2 sin Op^ d\pi I |p2 P sin 0p^d\p2 \ 



Pi 

9^1 



2-Er, 



./PDF 



9^2 



/flux(9zi , 922) /it 



xMti ^t2 ^-^Wi dMt^ dMt^ dMwi dMw^, 



(15) 



where the remaining integrations are over the invariant 
masses of the t quarks and the W bosons, and the mag- 
nitude of the h quark momenta. 

The cross section as a function of Mt is expressed as a 
six-dimensional integral; this integration is performed nu- 
merically using the VEGAS l39| algorithm as implemented 
in the GNU Scientific Librarv|40j. 



of the Zjj system is zero and instead integrate over the 
unknown pT which arises from additional softer jets and 
unclustered energy. We express this as rx and r^, compo- 
nents of the recoil in the x and y axes, respectively. We 
integrate over dr^dry using uncorrelated Gaussian priors 
in X and y with widths of 12 GeV/c^, fuTF{'fx,'>'y), as 
extracted from simulated samples. The differential cross 
section can be expressed as 



ddz 



33 



dx (2 



\ \ f ^ 



EliEl2Ep^Ep^ 2 
W{Ep,,E,,)W{Ep,,E,.,)\p,\^\p2\^ 



\qziqz2\\sin{(j)j. 



I^'J2, 



where pi , p2 are the four- momenta of the final state par- 
tons which lead to creation of extra jets, Li,L2 are the 
four-momenta of the final state leptons, and qi, 52 sue the 
four-momenta of incoming partons. 



VII. BACKGROUND LIKELIHOODS 



B. WW + 2 jets 



We calculate the per-event differential cross section for 
the three largest sources of background: Z/j* + 2 jet 
process {Zjj) where the Z decays directly to electrons or 
muons, the WW + 2 jet process (WWjj) and the -I- 3 
jet processes (Wjjj) where one jet is misidentified as a 
lepton. WZ with associated jets and Z tt with two 
jets have a small overall contribution to the sample and 
are not directly modeled. 

The major background processes can not be well de- 
scribed using a small number of diagrams. We therefore 
adapt routines from ALPGEN which make effective ap- 
proximations to evaluate the matrix elements for these 
processes. The ALPGEN routines are a function of the 
spin and color configurations of the initial and final state 
partons as well as their momenta. We employ a statis- 
tical sampling over the spin and color configurations to 
numerically evaluate the averaged M . 

The final measurement is calibrated using fully realistic 
Monte Carlo events, which will incorporate the effects of 
these approximations. 



A. Zf-y* + 2 jets 

We employ the set of assumptions as described in Sec- 
tion IV and use transfer functions as defined and derived 
in Section V to connect the parton-level quantities to 
observed quantities. We use ALPGEN's Z + 2p matrix 
element, which models the associated production of two 
light quark (w, d, c, s) and gluon jets. As the number of 
unmeasured quantities is fewer than in the case of the sig- 
nal likelihood, we can relax the assumption that the px 



The production of W pairs with associated jets is mod- 
eled in a similar fashion. We use ALPGEn's WW + 2p ma- 
trix element, which models the associated production of 
two light quark {u,d,c,s) and gluon jets. After making 
the assumptions listed in Section IV, we choose to trans- 
form the phase-space by introducing the invariant masses 
of the intermediate W bosons. We express all the par- 
ton variables except the neutrino momenta in spherical 
coordinates. Integrating over delta functions for lepton 
energy, jet angles, and total conservation of momentum 
gives 



da 



WWj] 



dx (2 



_\Ml_ 



a,b 

\pi\'^smep^d\pi\ |p2p sin 6*^2 d|p2 



2Ep-^ 



2Ep 

fa I Izi j ra 
"/PDF 7n ) JPl 



<1Z2 



Ei^2 ^L\^L2 \ -^beam J V -^becim 

X /fiux('72i , qz2W{Ep, , Ej, ) W{Ep^ , Ej^ ) 
X AfvKi ^W2 dMwi dMw^dvizdv2z, 

(17) 



where Li, L2 are the measured four-momenta , li, £2 are 
the parton-level four-momenta of the final state leptons, 
vi , V2 are the four-momenta of the final state neutrinos, 
Pi, p2 are the four-momenta of the final state partons 
that lead to creating extra jets, and gi, (72 are the four- 
momenta of the incoming partons. The sum runs over 
the incoming parton fiavors. 
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The final integration is performed over the the mo- 
menta of the partons which lead to jet production, the W 
boson invariant masses, and the z components of neutrino 
momenta. Transformation of the space requires solving 
a coupled system of equations to express the neutrino 
energies in terms of the W masses. 



C. Modeling of Backgrounds with Jets 
Misidentified as Leptons 

Events in which a jet is misidentified as a lepton can 
be modelled with the process pp Wjjj — > lyjjj- We 
use alpgen's + 3p matrix element, which models the 
production of three light quark (u, d, c, s) and gluon jets. 
Using this process as the basis for the model, we sum 
over the probability that either reconstructed lepton is a 
misidentified jet. Calculation of the cross-section follows 
the style of the other calculations above. 

Making the standard assumptions, and integrating 
over all delta functions gives 



Wjjj 



E d\pi\d\p2\d\p3\diy^ 

X W{Ep, , E,, )W{Ep, , Ej,)WiEp, , E,, ) 
X \pi\'^\p2f\p3fsin9j^ sin 6'j2 sin 6*^3 



ra I Izi \ rb ( 1Z2 
^ JPDF \ p I JmF I p 



\Ml 
2 

(18) 



where Pi,P2,P3 sltc the four- momenta of the final state 
partons which lead to creation of extra jets, L is the four- 
momenta of the final state lepton, v the four-momentum 
of the final state neutrino, and 51,(72 are the four- 
momenta of incoming partons. Jets which are misiden- 
tified as a lepton have the large majority of their mo- 
mentum carried by a single leading object (tt" or tt^); 
therefore, we assume that the jet identified as a lepton 
carries the momentum of the parton. 



VIII. TOP QUARK MASS RECONSTRUCTION 

In this section, we discuss the combination of the pcr- 
event differential cross section calculations for signal and 
background into a joint probability for a sample of events, 
the procedure for mass extraction, and calibration of the 
method using simulated events. 



A. Posterior Probability 

We express the individual event probability density in 
X as a sum of the signal and background probabilities 
with their respective fractions, see Eq. 4. The signal 
and background fractions are expressed in terms of the 



number of expected signal events Xs{Mt) and background 
Af, events: 



PsiMt) = 



A. (Aft) 
A,(Mt) + Af, 



(19) 



where A;, = J2i ■ The expected signal Xs{Mt) is calcu- 
lated relative to a reference point Mq and extrapolated 
to other masses using the mass dependence of the total 
accepted cross section as{Mt)es{Mt): 



XsiMt) = X,{mo) 



a,{Mt)e,{Mt) 
a,(Mo)e,(Mo)^ 



(20) 



where as{Mt) is the total production cross section p31j 
and €s{Mt) is the acceptance measured in Monte Carlo 
events. 

For an individual event x, P(x|A/t) is a likelihood in 
Mt ■ The posterior probability density in Mt is the prod- 
uct of a flat prior probability and the product of the 
individual event likelihoods. The measured mass. Aft, is 
chosen as the expectation value of the posterior probabil- 
ity to avoid potential fluctuations in the position of the 
maximum probability: 



Aft 



/ dAftAftP(x|Aft) 
/ dAftP(x|Aft) 



(21) 



The measured statistical uncertainty, AAft , is the stan- 
dard deviation of the posterior probability. 



AAft 



/ / dMtM^P{-x\Mt) 



J dMtP(MMt) 



[Mtf 



(22) 



B. Calibration 

We have searched for any potential biases on the ex- 
tracted mass or its uncertainty due to our fitting proce- 
dure. We parameterize this bias using a linear correction 
factor to the measured mass consisting of an offset Afo 
and a slope sj\/t , 

Aft = 178.0 GcV/c2 + {Mt - Afo)/sA/, , (23) 

and a simple scale factor, 

AAft^^AX AAft/5A/,, (24) 

to the measured statistical error in data, where 
Afo, sa/j , Sa are extracted from ensembles of Monte Carlo 
experiments. In the following sections, we study the cal- 
ibration of the method in Monte Carlo experiments with 
only t< events as well as in Monte Carlo experiments with 
both t< and background events. 
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FIG. 4: Top, mean of measured Mt in Monte Carlo exper- 
iments with only tt events. Bottom, fitted widths of pull 
distributions from the same Monte Carlo experiments. 



FIG. 5: Correlation between the mean pr of the tt system 
and the number of reconstructed jets in pythia simulated tt 
events. 



by a factor of Sa = 1-4, independent of top quark mass. 
This is due to the several simplifying assumptions made 
in the formulation of the probability expression. As a 
result, it does not contain a description of some effects 
which cause smearing of the extracted mass. The most 
important of these are the assumption that all jets come 
from b quarks, rather than from initial or final state QCD 
radiation, that jet angles accurately give the parton an- 
gles, and that lepton energies are perfectly measured. We 
use Monte Carlo experiments to measure the effects of 
these simplifying assumptions on the measured top quark 
mass. This is done by measuring our ability to extract 
the top quark mass in Monte Carlo experiments where 
these assumptions are violated to increasing degrees. 



C. Monte Carlo Experiments with tt Events 

We first construct Monte Carlo experiments of ti 
events with varying Mt, each generated by herwig [33 |. 
and a simulation of the CDF II detector. The Monte 
Carlo experiments contain a number of tt events drawn 
from a Poisson distribution whose mean corresponds to 
the number of tt events expected in a sample at the given 
mass. No background events are included in these exper- 
iments, and the pi,. in Eq. 4 are therefore identically zero. 
Figure 0] shows the mean measured Mt and the width of 
the pull distribution at each mass, where a pull is defined 
as {Mt — Mt)/ /S.Mt for each Monte Carlo experiment. 

The mean measured top quark mass in these Monte 
Carlo experiments shows no evidence of bias. A linear 
fit to the extracted mass yields Mg = 178.0 GeV/c^ and 
Sjvft = 1-00. The pull width, however, indicates that the 
extracted statistical error is consistently underestimated 



1. Jet-Parton Assignment 

In samples of simulated events which pass selection re- 
quirements, 70% of events contain two reconstructed jets 
whose axes lie within a cone of Ai? = ^/[Kqf^+{K(f)^ < 
0.7 of unique b quarks from the top quark decay. Monte 
Carlo experiments using this subset of events have a sig- 
nificantly smaller pull width, 1.20. This suggests that 
events in which the assumption of correspondence be- 
tween jets and 6 quarks is violated contribute significantly 
to non-unit pull widths. 

We note that the largest source of incorrectly assigned 
jets is initial state radiation, where a hard emission from 
the incoming quark or gluon gives rise to a jet. In sim- 
ulated events, a strong correlation is seen between the 
number of reconstructed jets and the pr of the tt sys- 
tem, which is a measure of the total recoil against initial 
state radiation; see Figure |S1 
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2. Lepton Resolution 



Though lepton energies are weU measured by CDF, 
electrons and muons are measured by different subdetec- 
tors. The energy of electrons at high Ex is very well 
measured by the calorimeter |14| : 



_ 13.5% 



2.0%. 



(25) 



The momentum of muons is measured by the central 
tracker, whose resolution begins to degrade at large 

PT M- 



=0.0011 •pT(GeV/c). 

PT 



(26) 



Monte Carlo experiments formed using events in which 
jets are matched well to b quarks and events containing 
only electrons have a pull width of 1.10. Similar Monte 
Carlo experiments using only muons have a pull width 
of 1.20, indicating that muon momentum resolution con- 
tributes to pull widths greater than unity. Electrons and 
muons which arise from W — > ri^r decays are not well de- 
scribed by the matrix-element so they contribute to the 
pull width as well. 



3. Jet Angle Resolution 

The jet angle resolution is finite, though it is sig- 
nificantly more precise in comparison to the jet en- 
ergy resolution. Figure El shows the angular distance, 
R = -y/ (Ar;)^ 4- (A(/))^, between reconstructed jets and 
the closest b quarks in fully simulated tt events passing 
the selection criteria in Section IIIII Jets which are not 
matched to either of the b quarks from top quark decay 
(no b quark within R < 0.7) are likely due to initial state 
radiation, as described above. 

To isolate the effect of the jet-angle resolution, we ex- 
amine a subset of events from a fully simulated sample 
with Mt = 178 GeV/c^. To remove the effects of lepton 
resolution and jet-parton matching as isolated above, we 
require well measured leptons (p^^*°" _ p-reconstructed ^ 2 
GeV/c) and matched jets (i? < 0.4). Events of this type 
have a negligible rate of jet-parton misassignment. Monte 
Carlo experiments with this subset of events have a pull 
width of 1.1, consistent with experiments described above 
which use only electrons and matched jets. Further tight- 
ening the R requirement reduces the pull width to 1.0. 

The effect on the pull width of jet angle resolution is 
confirmed by observation of a w 10% increase in pull 
width when jet angles are smeared in otherwise perfectly 
measured parton-level events. 
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FIG. 6: Angular distance, R, between a reconstructed jet and 
the closest b quark in tt events simulated with herwig. The 
width of the distribution demonstrates the angular resolution. 
20% of jet-parton pairs have R > 0.7 (not shown), coming 
from jets with no corresponding b quark. 



D. Monte Carlo experiments with tt events and 
background events 

To more fully describe the sample of events we expect 
in the data, we construct Monte Carlo experiments with 
ti events as well as events which model the background 
sources listed in Table Each Monte Carlo experiment 
contains numbers of events drawn from Poisson distribu- 
tions whose means are given by the expected number of 
events from each given source. 

To extract the mass, we use the full probability expres- 
sion of Eq.01 including terms which model the production 
of three background sources. Figure [T] shows the mean 
extracted mass and width of pull distributions in these 
Monte Carlo experiments. A linear fit to the extracted 
mass yields 



Mo = 177.2 ± 0.21 GeV/c^, sm, = 0.84 ± 0.02 

as defined in Eq. 1231 These parameters indicate a small 
bias due to lack of modeling of the minor backgrounds 
(WZ and Z tt) and imperfect descriptions of major 
backgrounds. If the background probability calculations 
are not included, Mq becomes 175.8 and SMt decreases to 
0.72. This demonstrates that the inclusion of the back- 
ground probabilities does improve the sensitivity of the 
method. The width of the pull distributions in fully real- 
istic Monte Carlo experiments using the full probability 
expression is fit to a straight line to extract the error 
scale factor, 



Sa = 1.51 ±0.02. 
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IX. SYSTEMATIC UNCERTAINTIES 
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FIG. 7: Top, mean of measured Mt in Monte Carlo experi- 
ments with ti events as well as background events. Bottom, 
fitted widths of pull distributions from the same Monte Carlo 
experiments. 



This scale factor is larger than the scale factor of 5a = 
1.4 measured in experiments with no background events. 
This is due to the presence of unmodeled background 
events {WZ and Z tt), which smear the extracted 
mass without inflating the measured statistical error. 

Uncertainty on the measured scale factor would con- 
tribute a systematic error to the statistical error of the 
measurement. The measured scale factor is stable to 
within 5% under variations of the sources of system- 
atic uncertainty in the modeling of the calibration Monte 
Carlo described in Section IX below. 

After application of mass and error corrections, the 
final measured mass is unbiased and the pull distribution 
is well described by a Gaussian; sec Figures |S1 and O 



The measurement is calibrated using Monte Carlo sim- 
ulated events. Therefore, the majority of systematic un- 
certainties come from uncertainties in this modeling of 
the data. In this Section, we describe the facets of the 
simulation which may not accurately describe the ob- 
served data and estimate the effect on the measurement. 

To measure the size of the impact of each uncertainty, 
we perform Monte Carlo experiments using a pool of sim- 
ulated events in which a feature of the events has been 
modified. By extracting the average measured mass, we 
can determine the typical shift due to these features. 



Systematic Uncertainties Due to Jet Energy 
Scale 



The largest source of systematic uncertainty arises 
from potential mismodcling of the jet energy measure- 
ment, through uncertainties in the various corrections 
applied to the measured jet energy These jet en- 

ergy corrections involve knowledge of the absolute energy 
scale, energy loss outside the jet search cone AR, the 
non-uniformity in response of the calorimeter as a func- 
tion of 77, effects from multiple pp collisions and energy 
deposition from the underlying pp event. These uncer- 
tainties are shown as a function of jet Pt in Figure [TUl 
A systematic uncertainty is estimated for each jet en- 
ergy correction by performing Monte Carlo experiments 
drawn from simulated signal and background events with 
±1 standard deviation in correction uncertainty, and tak- 
ing the half-difference in mean reconstructed top quark 
mass between the two results. The uncertainties from 
each energy correction are then added in quadrature to 
arrive at a total systematic uncertainty on the jet energy 
scale. No strong dependence was observed as a function 
of the top quark mass. 

Since the above jet energy corrections are developed 
from studies of samples dominated by light quark and 
gluon jets, additional uncertainty occurs from extrapolat- 
ing this procedure to b quarks. The resulting systematic 
effect on jet energy is considered to stem from three main 
sources: uncertainty in the b jet fragmentation model, 
differences in the energy response due to semi-leptonic 
decays of b hadrons, and uncertainty in the color flow 
within top quark production and decay to b jets j4lj . 
As in the jet energy scale uncertainty, Monte Carlo ex- 
periments arc performed on events where the b jet ener- 
gies have been altered by ±1 standard deviation for each 
uncertainty, and the resulting half-differences added in 
quadrature to estimate the total systematic uncertainty 
due to b jet energy uncertainty. 

Table Unl lists the uncertainties due to each correction. 
The sum in quadrature of the uncertainties in Tabic IIIII 
yields a flnal jet energy scale systematic uncertainty of 
2.6 GcV/c2. 
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FIG. 8: Left, distribution of measured Mt in Monte Carlo experiments with ti events at Mt = 165 GeV /c? and background 
events. Measured statistical error, center, and pull distribution, right, are also shown. All mass and error corrections have been 
applied. 




FIG. 9; Left, distribution of measured Mt in Monte Carlo experiments with ti events at Mt = 178 GeV /c? and background 
events. Measured statistical error, center, and pull distribution, right, are also shown. All mass and error corrections have been 
applied. 



TABLE III: Systematic uncertainties on the top quark mass 
due to each jet energy systematic uncertainty. 

Jet Energy Systematic 5Mt (GeV/c^) 



Modeling of hadron jets (absolute scale) 


1.1 


Modeling of parton showers (out-of-cone) 


2.2 


Response relative to central calorimeter 


0.8 


Underlying event and multiple interactions 


0.1 


Modeling of & jets 


0.5 


Total jet energy systematic uncertainty 


2.6 



B. Systematic Uncertainties Due to Backgrounds 



There are two major sources of systematic uncertainty 
that result from the presence of background. One is the 
limited number of events to model the background and 
calibrate Monte Carlo experiments. The other is the un- 
certainty in modeling major sources of background. 



1. Background Statistics 

We estimate the sensivity of the measurement to the 
limited number of events available to model the back- 
ground processes in the calibrating Monte Carlo experi- 
ments. We perform Monte Carlo experiments using dis- 
joint subsets of statistically independent event samples 
for each background process. The root mean square of 
the difference between the mass measured in each pair 
of experiments is an estimate of this uncertainty and is 
measured to be 1.2 GeV/c^, driven largely by the limited 
amount of data available for modeling events with a jet 
misidentified as a lepton. 



2. Background Modeling 

We estimate the sensitivity to the modeling of the two 
largest backgrounds, Drell-Yan and misidentified leptons. 

The contribution to the sample from Drell-Yan produc- 
tion comes from events in which there is an apparently 
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FIG. 10: Systematic uncertainties for several sources in the jet energy scale, as a function of the corrected pr of jets {p'r'^^) in 
the central (0.2 < \ri\ < 0.6) calorimeter. 



large missing transverse energy due to mismeasurement. 
To gauge the sensitivity to events on the tail of the dis- 
tribution, we vary the composition of the Monte Carlo 
experiments by either enhancing events or suppressing 
events on the tail. We assign weights to the events pro- 
portional to the measured missing transverse energy, and 
to its inverse: 



W+ (X ^T, Wo = 1, W- OC 

This prescription conservatively describes the uncer- 
tainty in modeling of the missing transverse energy in 
events from DrcU-Yan production, sec Figure ^2 To 
determine the systematic uncertainty arising from the 
shape of the Drell-Yan background, we compare the mass 
resulting from Monte Carlo experiments using and 
W- weights to that obtained using unit weight (wq). The 
larger of the two differences is taken to be the systematic 
uncertainty and is measured to be 0.4 GcV/c^. 

The background resulting from events in which a jet is 
misidentificd as a lepton is very difficult to describe accu- 
rately with simulation, as it is very sensitive to the details 
of the detector performance. To avoid issues of model- 
ing, the events which imitate this background are drawn 
from the data themselves. These events are selected with 
the same requirements as for the candidates, except for 
a looser requirement on one of the leptons, and then are 
weighted by the probability that the loose lepton would 
pass lepton identification requirements. These weights 



are measured as a function of the pr and isolation of the 
fake candidate using a sample dominated by jets; each 
candidate has its own weight (w) and uncertainty (Aw). 
To gauge the sensitivity to the calculation of the fake 
rates, we vary the fake rates in two ways. First (a), we 
enhance those events with fake probability greater than 
the mean (w > w) to exaggerate their effect; second (6), 
we enhance events with fake probability smaller than the 
mean(u' <w), to exaggerate their effect: 

Mode w > W w <w 

(a) w —^ w + Aw w w — Aw 

(b) w —* w ~ Aw w w + Aw 

We take the difference between the mass obtained us- 
ing (a) and that obtained using (b) as the systematic un- 
certainty, 0.6 GeV/c^, and the total shape uncertainty to 
be the two (Drell-Yan and fakes) summed in quadrature, 
0.8 GeV/c2. 

C. Parton Distribution Function Uncertainties 

To propagate uncertainties in the parton distribution 
functions (PDFs) to the mass measurement, we reweight 
a sample of simulated PYTHIA events according to 20 sets 
of ±1ct uncertainty eigenvectors based on CTEQ6M [i^ . 
The uncertainty using each of the eigenvectors is added 
according to the standard prescription to yield the total 
uncertainty of Apjjp = 0.55 GcV/c^. 
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Missing Transverse Energy [ GeV ] 



FIG. 11: Comparison of the missing transverse energy in 
dilepton events in data, pythia simulated Zj-^* — > II events 
and simulated events weighted as described in the text to en- 
hance or suppress events with large missing transverse energy. 
Top shows [III, data and bottom shows ee data. 



The two samples MRST72 and MRST75 use differ- 
ent values of Aqcu (228 and 300 MeV) to calculate the 
PDFs. The difference in the extracted mass from Monte 
Carlo experiments using these two samples provides a 
measure of the sensitivity to the uncertainty in A^qcd- 
This uncertainty is measured to be 0.85 GeV/c^. 

The leading order pythia and herwig generators con- 
tain ^5% gg events in their initial state. In pp collisions 
at ^/s = 1.96 TeV , approximately 15% of the ti pairs are 
produced from gg annihilation. The matrix element used 
in the likelihood calculation describes ti production from 
qq annihilation only. To measure the sensitivity to the 
initial state, we vary the fraction of gg initial states in 
the Monte Carlo experiments from 0% to 15%. We take 
the difference of the measured mass using 0% and 15% 



gg events as a systematic uncertainty, which is measured 
to be 0.35 GeV/c2. 

The total PDF uncertainty is obtained by adding 
the above three uncertainties in quadrature, and is 
1.1 GeV/c2. 



D. Systematic Uncertainties Due to Initial and 
Final State Radiation 

The calculation of the signal probability does not con- 
tain a description of initial (ISR) or final state radiation 
(FSR), which may contribute significantly to the fraction 
of mismeasurcd events. The rate of additional jets from 
initial state radiation can be estimated in Z+jets events, 
as there are no significant contributions from final state 
radiation; the radiation is found to depend smoothly on 
the Drell-Yan mass squared ^j, and can be examined 
over a broad range of energies, extending up to the range 
of ti production. To measure the uncertainty due to im- 
perfect knowledge of the rate of radiation, we examine 
the measured top quark mass in samples where the Monte 
Carlo generator parameters are varied by very conserva- 
tive amounts based on the studies in the Z-|-jets events. 
We measure the uncertainty due to ISR to be 0.5 GeV/c^. 

Final state radiation can be probed in the same man- 
ner, as it is described by the same showering algorithm. 
We find the uncertainty due to FSR to be 0.5 GeV/c^. 



E. Other Systematic Uncertainties 

To account for a possible bias of the Monte Carlo mod- 
eling, we measure the difference in top quark mass as ex- 
tracted from HERWIG and PYTHIA samples. This amounts 
to 0.8 GeV/c^. 

In addition, uncertainties in the fitted response cali- 
bration described in Section IVllI Bl must be taken into 
account. We measure this by varying the response by Itr 
within the statistical uncertainties of the fit and measure 
the difference in extracted top quark mass. We find this 
uncertainty to be 0.4 GeV/c'^. 

The response calibration is derived from Monte Carlo 
experiments in which the expected number of ti events is 
calculated using the theoretical cross section as a func- 
tion of Mt- The 10% uncertainty on the cross section is 
propagated to the final mass measurement, yielding an 
uncertainty related to the actual sample composition of 
0.3 GeV/c2. 

Finally, we verified that the measurement is not sen- 
sitive to effects beyond leading order by measuring 
the mass in Monte Carlo experiments constructed with 
events with = 175 GeV/c^ generated by MC@NLO 
li^ l which includes next-to-leading order effects; the mean 
extracted mass in these experiments was 175.8 ± 0.8 
GeV/c2. 
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TABLE IV: Summary of systematic uncertainties. 



Source 5Mt (GeV/c^) 

Jet Energy Scale 2.6 

Background statistics 1.2 

Background modeling 0.8 

PDFs 1.1 

FSR modeling 0.5 

ISR modeling 0.5 

Generator 0.8 

Method 0.4 

Sample composition uncertainty 0.3 
Tbtal 



F. Summary of Systematic Uncertainties 

The systematic uncertainties measured in the previous 
sections are summarized in Table Hvl The total system- 
atic uncertainty is 3.4 GeV/c^ of which the single largest 
source is the jet energy scale, contributing 2.6 GeV/c^. 

X. MEASUREMENT 

The data sample in / Ldt = 340 pb~^ contains 33 can- 
didate dilcpton events. These candidates have individual 
likelihoods as seen in Figure E| From the joint proba- 
bility, we extract the uncorrected, unsealed mass, 

Mt = 166.4 ±3.4(stat.) GeV/c^ 

After applying corrections to the central value and sta- 
tistical uncertainty as derived in Section VIIID, the final 
result is 

Mt = 165.2 ±6.1(stat.) GeV/c^ 

The corrected joint probability curve can be seen in 
Figure [T^ In Monte Carlo experiments where Mt = 165 
GeV/c^, 17% of the uncertainties are smaller than this 
value, see Figure PPT] 

XI. CONCLUSION 

We report the first application of a matrix-element 
based method to the measurement of the top quark mass 
in ti events containing two leptons. We measure 

Mt = 165.2 ± 6.1(stat.) ± 3.4(syst.) GeV/c^ 

which is the most precise determination to date of the top 
quark mass in dilcpton events. This result is consistent 
with recent measurements of the mass in this channel at 
CDF using template methods. 
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FIG. 12: Likelihood in top quark mass for the 33 candidate 
events in 340 pb^^ of run II CDF data. Horizontal axis is 
top quark mass over the range of Mt = 130 to 220 GeV/c^. 
Vertical axis is shown on a linear scale between and 1. 
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with measurements in run I from CDF, 
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FIG. 13; Joint probability density for the 33 event data sam- 
ple in 340 pb"^ of run II CDF data as a function of the top 
quark mass. The bias and error corrections have been applied. 
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FIG. 14: Distribution of expected uncertainties for Mt = 165 
GeV/c^ in Monte Carlo experiments. The measured uncer- 
tainty is shown as the line; 17% of Monte Carlo experiments 
yield a smaller uncertainty. 



170.1 ±6.0(stat.)±4.1(syst.) GeV/c^ |4£ 



Mt = 167.4 ± 10.3(stat.) ±4.8(syst.) GeV/c^ @, 
and D0, 

Mt = 168.4 ± 12. 3(stat.) ± 3.6(syst.) GeV/c^ 0. 

This measured value is smaller than the current preci- 
sion measurements of the mass in single lepton events, 

Mt = 173.5l|^(stat.) ± 1.3(syst.) GeV/c^ |4|. 

A global combination of the most precise measure- 
ments [4^ . however, suggests that current discrepancies 
are consistent with statistical fluctuations. 

Extrapolating the use of the method to a future top 
quark dilepton data sample from the Tevatron corre- 
sponding to an integrated luminosity of J Ldt = 4 fb~^, 
the expected statistical uncertainty of this technique in 
Monte Carlo experiments is 2.5 GeV/c^, for Mt ~ 178 
GeV/c^. In this regime, uncertainty in the jet energy 
scale would be the dominant source of uncertainty. 
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